Self-organization of irregular NEM vibrations in multi-mode shuttle structures 
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We investigate theoretically multi-mode electromechanical "shuttle" instabilities in DC voltage- 
biased nanoelectromechanical single-electron tunneling (NEM-SET) devices. We show that initially 
irregular (quasi-periodic) oscillations, that occur as a result of the simultaneous self-excitation of 
several mechanical modes with incommensurable frequencies, self-organize into periodic oscillations 
with a frequency corresponding to the eigenfrequency of one of the unstable modes. This effect 
demonstrates that a local probe can selectively excite global vibrations of extended objects. 
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In a nanoelectromechanical single-electron tunneling 
(NEM-SET) device mechanical vibrations and single- 
electron tunneling events are coupled on the nanome- 
ter length scale. Ten years ago Gorelik et al. [if sug- 
gested that such a coupling of mechanical and electrical 
degrees of freedom could lead to a "shuttle instability" 
and to a novel shuttle mechanism for charge transport 
through a DC voltage-biased NEM-SET. The theory tU] 
was developed for a rigid nanometer-size metal cluster 
suspended by elastic links between a source and a drain 
electrode so that only one mechanical degree of freedom 
was involved. Internal cluster vibrations were assumed 
to have much higher frequencies than the relatively low- 
frequency center-of-mass vibration and were therefore ig- 
nored. In the proposed shuttle regime charge transport 
between source and drain is mechanically assisted by pro- 
nounced center-of-mass vibrations of the cluster (for re- 
views see 0, 

In this Letter — motivated by the recent experimental 
discovery of the shuttle instability predicted in Ref. [ij 
and of the related self excitation of radio-frequency me- 
chanical vibrations of gold-capped silicon nano-pillars Q 
— we generalize the theory of the shuttle phenomenon 
to include extended structures with (many) internal me- 
chanical degrees of freedom. Our objective is to stim- 
ulate further experimental work on NEM-SET systems 
involving, e.g., suspended carbon nanotubes and other 
extended molecules. In many such systems internal me- 
chanical degrees of freedom cannot be ignored and there- 
fore a theory of shuttling in the multi-mode regime needs 
to be developed. This will be done in what follows. 

Our main result is that when self excitation involves 
more than one mode, strong mode-mode interactions 
caused by the non-linear coupling of the mechanical 
vibrations to the electronic subsystem results in self- 
organization of the vibrations. This implies that ini- 
tially quasi-periodic vibrations are transformed into pro- 
nounced periodic vibrations as the shuttle instability de- 
velops towards a steady state. The frequency of the 
steady-state vibrations depends crucially on initial con- 
ditions and corresponds to the eigenfrequency of one of 
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Figure 1: Sketch of the model system considered. The elec- 
tron tunneling rate rL(a;) between the STM-tip and the center 
of a suspended carbon nanotube depends on the tube posi- 
tion, while the tunneling rate Fr between tube and electrodes 
is constant. Tunneling may lead to self-excitation of multi- 
ple nanotube bending modes, which self-organize into perio- 
dic, single-mode oscillations (the vibration amplitude shown 
is exaggerated). 



the unstable modes. 

We will consider the particular realization of a multi- 
mode shuttle structure sketched in Fig. [T] The sketch 
shows a NEM-SET device based on a carbon nanotube 
(CNT) suspended over a trench so that a segment of 
the tube is free to move in response to external forces. 
Charge is injected into the suspended nanotube through 
the tip of a scanning tunneling microscope (STM) posi- 
tioned above the center of the tube Q as in the exper- 
iment of LeRoy et al. Q. For this very system a cur- 
rent driven nanoelectromechanical instability has been 
predicted and theoretically shown to result in the onset 
of pronounced CNT bending mode vibrations involving 
one or more modes 0, 0|- Below we will consider the 
detailed development of the predicted instability when 
several bending modes are unstable in order to find the 
resulting steady state. While a nanoelectromechanical 
instability involving these modes has not yet been ob- 
served experimentally, bending mode vibrations of sus- 
pended CNTs have been detected in devices of the type 
studied here . We note in passing that the elec- 

tromechanical coupling to other types of CNT vibration 
modes have also been studied (gI. [ill [l^, [isl. [l^ . 

We use continuum elasticity theory to describe the 
nanotube mechanics 15 , 3] ■ The set of normal modes for 
the bending mode vibrations is characterized by the fre- 
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^^'^ , where E is Young's modulus of 



quencies =- jjsy 
the CNT, Iz is the area moment of inertia of the cross sec- 
tion, p is the mass density, S is the cross section area and 
the coefHcients c„ = 22.4, 61.7, 120.9, 199.9, 298.6, ... are 
obtained by solving the equation cos cosh ^/c^ = 1 
17| . The displacement x{t) of the center of the CNT 
(see Fig. [1]) can be expressed as a sum of normal-mode 
amplitudes x„{t), where only modes with n = 1, 3, 5, 
. . . contribute due to symmetry. Other modes are inert 
and can be neglected. Each normal mode coordinate Xn 
obeys Newton's equation in the form 



Xji -\- ^Xji -j- UJj^Xrt 



q£/M, 



(1) 



where the force term depends on the average net charge 
q on the CNT and the effective electric field £ that acts 
on the tube as a result of the applied bias voltage 
The damping coefficient 7 in Eq. ([T]) is a constant in 
our model and therefore affects all modes in the same 
way. It is related to the quality factor of the vibration 
modes through the relation Qn = tOn/l- Experimentally 
Q factors of order 1000 have been reported for nanoscale 
CNT based resonators in vacuum, [18|, while Q factors 
as high as 10*— 10^ have been found for somewhat larger 
Si dll and SiC resonators 

An important parameter for the system under consid- 
eration is the ratio lu/T between the typical CNT vibra- 
tion frequency w and the characteristic rate of electron 
tunneling T. Having in mind that the mechanical vi- 
bration frequency is extremely low on an electronic scale 
{u) ~ 10^ — 10^ s""'^ for L ~ 1 /im) we consider here the 
case Lo <^T. In addition we focus on the single-electron 
shuttling regime, where only one electron can be accumu- 
lated on the CNT due to Coulomb blockade of tunneling. 
Introducing the probability p to find an extra electron on 
the nanotube the average excess charge is g = ep while 
the kinetic equation for the time evolution of the proba- 
bility p can be written as 



p = -T{x)p + Ti,{x). 



(2) 



Here T{x) — Yi^{x) + Fr with x — S,ia;„ while FL(a;) — 
Fq exp(— x/A) is the rate of electron tunneling across the 
STM-CNT junction, so that the typical rate of electron 
tunneling is F = Fq + Fj^. The characteristic length A is 
known as the tunneling length. 

Equations 111 and ([2]) describe the coupled nanoelec- 
tromechanical dynamics of the CNT-based NEM-SET 
device. In order to stay in the single-electron shuttling 
regime one has to apply a small enough bias voltage: 0.1 
V is a typical value for a CNT of length L ~ 1 /im. There- 
fore, it is reasonable to consider the limit of weak nano- 
electromechanical coupling where the tunneling rates are 
only weakly modified by the electrostatic force induced 
by an excess charge q = e. This condition holds for 
a small enough voltage bias V (and a correspondingly 



small effective electric field £) when the shift of the nan- 
otube equilibrium position due to a single excess charge, 
dn — e£/{Muj^), is small on the scale of the tunneling 
length, i.e. when (i„/A <C 1. In this limit the onset of 
a shuttle instability occurs independently in the differ- 
ent vibration modes Q and if cj ^ F the instabilities 



are soft [2l|, i.e. the amplitude An of stationary shuttle 



vibrations goes to zero as £ approaches a critical field 
£c = ATX'^M/e from above. For small enough positive 
values of £ — f c the vibration amplitudes An will there- 
fore be small compared to the tunneling length A and 
one may expand in the small parameter A„/A. By keep- 
ing third order terms one captures the most important 
non-linear effects. 

Now we analyze Eqs. ([T|) and Q in the limits discussed 
above. A formal solution of Eq. ^ is 



pW = E(-r(x)"'^0'"rL(:c)/F(a:) 



(3) 



m— 



which is a series expansion in the small parameter lo/T. 
If aj„ ^ F for all modes n it is sufhcient to retain 
only the first order term and substitute the truncated 
solution into Eq. ([1]). After expanding to third order 
in the displacements Xn, the resulting non-linear equa- 
tions for Xn are solved by choosing the Ansatz: Xn{t) — 
XAn{t) sm{ujnt -\-J{ja{t)) and by then averaging over the 
fast oscillations [2^. The remaining equations describe 
the slow time variation of the amplitudes A„ and phases 
Xn of the vibrations (A„, <ti w„). Since the oscillator 
frequencies w„ are incommensurable (see above) the am- 
plitude and phase equations are decoupled in the limit 
considered here. Then the amplitude equations are 



A-n — 0;nAn{Sn — A'^ 



(4) 



with 

Sn = 16 
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One notes from the expression of d„ that the product 
ujndn does not depend on n and that therefore Sn and 
in Eq. ^ do not depend on n to leading order in Wn/F. 
This will be used in the following analysis and allows us 
to write a„ = a and Sn — S. The relevance of corrections 
to these values will be discussed later. 

A complete analysis of how the solutions to Eq. (j4]) 
evolve is possible when two modes n, m are unstable. 
In this case one has two non-linear first order differen- 
tial equations and the stationary points can be classi- 
fied using standard techniques. The stationary points 
found and their classification are listed in Table [T] and 
the phase space is shown in Fig. [2l We find that the sta- 
tionary solution corresponding to finite amplitude vibra- 
tions of both modes is a saddle point and the solutions 
describing finite amplitude vibrations of one mode and 
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{^1,^3} 


Type 


{0,0} 




Repulsive fixed point 






Attractive fixed point 


{V5,0} 




Attractive fixed point 






Saddle point 



Table I: Classification of stationary points when two modes 
{n—1 and 3) are included and both are unstable {S > 0). Four 
stationary points where Ai,A3 > are found. Two of them 
are attractors and correspond to only one unstable mode, the 
other being stable. Which mode is unstable is determined by 
the initial conditions. 




Figure 2: (Color online) Stationary points when two modes 
(n and m) are unstable. Two attractors, indicated by (o), 
correspond to a finite amplitude of one mode while the other 
mode is suppressed. The stationary point marked with (x) 
is a repellor and the point indicated by (*) is a saddle point. 
The thick lines are separatrices that trajectories cannot cross. 
The separatrix An = Am ensures that if A„{0) > Am{0), this 
inequality hold for all times t. 



zero amplitude of the other are attractive fixed points. 
This means that depending on the initial conditions, one 
mode is selected to vibrate with finite amplitude while 
vibrations in the other mode are suppressed. Although 
a complete analysis of the general non-linear problem 
can not be done, a general statement for a system with 
N > 2 modes can be formulated for small 6. In this 
limit one finds that the N coupled non-linear equations 
^ are characterized by N stable fixed points in the 
vector space A — {An}. The evolution of the system 
in this space is represented by the motion of the point 
A{t) along a certain trajectory. The stable points are 
given by Ai{oo) = (A;(oo) = VS, An^i{oo) = 0), cor- 
responding to one mode vibrating with amplitude \/^ 
and all other modes having zero amplitude. More pre- 
cisely we have proven the following theorem: if the initial 
conditions are such that Ai(0) > An^i{0), then the am- 
plitude Ai{t) increases monotonically towards the final 
value Ai{oo) = ^/S while the other amplitudes An^i{t) 
decay exponentially towards zero as t oo. The asymp- 
totic estimate An^i{t) < VS exp{-at{Af{0) - ^^^(0))} 
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Figure 3: (Color online) Nanotube vibration amplitude as 
a function of time when three modes (n = 1,3,5) are un- 
stable {5 > 0) and included in a full numerical solution of 
Eqs. (HJ-Q with 7/1.^1 = Q-^ = 6.1 x 10"®, di/A = 0.01, 
wi/r = 0.0025 and initial conditions Ai{0) = 0.02, ^3(0) = 
0.01, ^5(0) = 0.03. These parameters were chosen to allow 
a comparison with the approximate result of Eq. Q which 
is plotted as dashed curves. The lower left panel shows the 
quasi-periodic oscillation of the nanotube center position just 
after the onset of the instability, while the lower right panel 
shows the regular vibrations that appear after all but the n=5 
mode amplitudes have been suppressed (see text). 



gives an upper limit to the decaying amplitudes. The 
theorem guarantees a remarkable feature of the shuttle 
instability in that it allows a selective amplification of 
one of the normal modes of the shuttle vibrations. Ini- 
tially quasi-periodic vibrations (for incommensurable fre- 
quencies ujn) are forced to self-organize into regular large- 
amplitude vibrations with a frequency corresponding to 
that of a single normal mode. 

We now resort to a numerical analysis in order to verify 
and generalize the analytical results obtained above. We 
begin by using the same parameter range as before, i.e. 
we choose parameters in order to allow a comparison be- 
tween the numerical and analytical results rather than to 
model realistic experiments. We will then show that the 
same type of behaviour obtains for more realistic sets of 
parameters, for which no analytic solutions is available. 
Figure [3] shows numerical and analytical results for the 
time evolution of the vibration amplitude when there are 
three different unstable mechanical modes {n=l, 3, 5). 
It is clear that even though all three modes initially in- 
crease their amplitudes, a single mode is selected and the 
other modes are suppressed at large times. Here the n=5 
mode initially had the largest deviation and hence the vi- 
bration frequency in the stationary state is given by lu^. 
Similar figures when modes n=l or 3 ends up with a fi- 
nite amplitude can be obtained by changing the initial 
conditions. 

So far we have only discussed cases when the final vi- 
bration amplitude is smaller than the tunneling length A. 
Numerically, we have also investigated the opposite situ- 
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Figure 4: (Color online) Time evolution of the vibration am- 
plitude obtained from a numerical solution of Eqs. ([T|l-(f2)) 
when three modes (n = 1,3,5) are unstable and 7/ciJi — 
= 5.5 X 10^'*,di/A = 0.7 and loi/F = 0.005. The 
initial conditions were ^i(O) = 0.005, ^3(0) = 0.0025 and 
^5(0) — 0.01. The large amplitudes make an approximate 
analysis based on Eq. Q invalid, but the the phenomenon of 
a selective excitation persists even in this case — the mode 
with the largest initial deviation is selected to end up with a 
large vibration amplitude. 



ation and found results that are qualitatively similar. As 
an example, a full numerical solution of Eqs. ([T]) and ([2]) 
for a case when the final amplitude exceeds A is shown 
in Fig. H 

Finally we discuss the coefficients a„ and (5„. In the 
previous analysis we used the approximation that these 
coefficients are the same for all modes, which is only valid 
when F 3> max u;„ and the dissipation mechanism affects 
all the modes in the same way. In general if one takes 
corrections of order w^/F^ into account or considers a 
mode-dependent dissipative mechanism, one has (5„ 7^ 
5m- This breaks the symmetry between the modes and 
leads to a situation where the final vibration state is not 
simply determined by the largest initial vibration mode 
amplitude. In this case the number of stable attractors 
can be smaller than the number of unstable modes. This 
implies that in some cases, even though a mode is initially 
unstable, it cannot reach finite amplitude vibrations since 
it is suppressed by another mode. 

In conclusion we have investigated electromechanical 
instabilities of different mechanical vibration modes of a 
suspended carbon nanotube. We have shown that the 
excitation mechanism considered leads to a selective ex- 
citation of a specific mode depending on the initial condi- 



tions. This demonstrates a way of using local tunneling 
injection of charge to probe the mechanics of extended 
nano-objects. The analysis presented here can be gen- 
eralized to apply to other multi-mode shuttle structures. 
The phenomenon of a selective excitation of a specific 
mode is general in the sense that even though several 
modes are unstable, only a single mode reaches a steady 
state with a finite amplitude. 
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